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In this paper, we study the spin-spin interaction between two electrons bounded in a quantum 
dot. The result shows that spin-spin interaction will cause a pair of spins precessing synchronously. 
If the two spins are parallel at initial time, the total spin oscillates as a cosine function. If the two 
spins are antiparallel at initial time, the total spin keeps zero. The precessing period is proportional 
to the cube of quantum dot radius. For a 2D round quantum dot with radius R — 50nm, the 
precessing period equals nearly 12/is. 

PACS numbers: 72.25.Rb , 03.65.Ge 

I. INTRODUCTION 

The spin dynamics in quantum dots attracts much attention in recent years, because the experiments indicate long 
spin relaxation time []) Q 9 m some kinds of quantum dots. This makes the quantum dot the promising candidate 
for qubit, which is the foundation of quantum computer. Many plans have been proposed for representing qubit and 
building unitary gates in quantum dots system. Some authors have researched the spin filling in a quantum dot [ij 
and proposed the qubit encoded by more than one electron spin, such as encoding a qubit in singlet and triplet states 
for two electrons in one dot or in two coupled dots [l| [f| . Experiments have shown the singlet-triplet transition can 
be electrically controlled QGUl- 

Furthermore, there are proposals for encoding a qubit in two specific spin states of 
three electrons in single dot [9J. In order to get a good representation of qubit and the corresponding gate sequences, 
it is important to understand the interaction on spins in quantum dots. Hyperfine interaction between the nucleus 
and electrons has been considered as the most important interaction on spins in some kinds of quantum dots (Toj [ll[ . 
It has been investigated by many authors [12] 13][14][l5j. 

The Coulomb interaction between bounded electrons has also been intensively studied. But according to quantum 
electrodynamics (QED), Coulomb interaction is only the non-relativistic approximation of electron-electron inter- 
action. Considering the relative nature of electron movement, spin-spin interaction, which can be regarded as the 
interaction between two magnetons, should be included. There is great difference between Coulomb interaction and 
spin-spin interaction, because the latter can exert forces on spins. Considering spin-spin interaction, only the total 
angular momentum conserves. Spin-spin interaction has been researched in condensed matter physics after the QED 
theory was discovered. Overhauser [161 ] in his early articles has discussed the spin relaxation caused by spin-spin 
interaction between Bloch electrons. He found that the effect was much smaller than that caused by spin-orbit inter- 
action and could be ignored. But in a quantum dot, the spin-orbit coupling is efficiently suppressed, while spin-spin 
interaction is enhanced because electrons are much closer to each other. On the other hand, the experiment about 
spin filling shows that there exist unpaired electrons in a multi-electron quantum dot. Then there exists the spin-spin 
scattering process caused by spin-spin interaction. 

In quantum dots, spin-spin interaction does not cause spin relaxation or dephasing. But it has non-trivial influence 
to the evolution of spins in quantum dots. Therefore it is obliged to study the effects of spin-spin interaction. 

The paper is organized as follows: in Sec. II we give the Hamiltonian of spin-spin interaction from QED theory 
Then we discuss the selection rule of this interaction in Sec. III. In Sec. IV and Sec. V, we study the evolution 
process of two parallel and antiparallel spins respectively. Sec. VI is a brief conclusion. 

II. SPIN-SPIN INTERACTION 

The Hamiltonian of QED theory can be expressed as 

H = Hj d3 P w P c L c P.s + J d 3 k{-gxx\k\)al x a KX -e J dPrftytyAp. (1) 
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Here we have omitted the terms creating or annihilating positrons. The Fourier decompositions of the field operators 
are 



/i7 3 k 
-,_= 5> M e*(k, A)e ikr + 4, A e"(k, A) e - tkr ). (2) 

By using the canonical transformation and noticing the relation g\\e^(\i, A) e w (k, A) = g^j,, one can get an effective 
electron-electron interaction 

/" e 2 

^/Z = H J ^ 3 P^ 3 P^ 3 k 2(2 ^ )3|k|2 u(p + k, s) 7 "u(p, s')u(p' - k, fl"')7/i«(p'. s > P +k, s c P,s' c p'-k, s '" c P',s'' ■ (3) 



S.S' , 8" .S 



Here we have used the low energy approximation uj p ps m. There are two terms in this effective Hamiltonian which 
exert forces on the spins. The first is the spin-spin interaction which is the coupling between one electron spin and the 
magnetic field caused by another electron spin. And the second is the spin-current interaction which is the coupling 
between the magnetic moment of one electron and the magnetic field produced by the translational motion of another. 
The spin-current interaction between bounded electrons can be neglected. The spin-spin interaction can be divided 
into two terms 

V s - S = V D + V ex . (4) 

The first term Vn represents flip of a pair of parallel spins, and can be called parallel interaction. The second term 
V ex represents flip of a pair of antiparallel spins, and can be called antiparallel interaction. According to Eq. ([3]), the 
parallel interaction can be expressed as 

v D = J d^ d w fjf;j^ ]2 cUA>-^'^ + (5) 

where | and j denote the spin-up and spin-down states in z direction respectively. And c Pi -[ denotes the annihilation 
operator of spin-up electrons. The antiparallel interaction can be expressed as 



f . , . ,-e 2 h 2 (\k\ 2 + k 2 z ) t + 
V ex = J d kd pd p 4m2c 2 eo(27r)3 | k |2 c P+k,T c p'-k,i c P'.T c P4- ( 6 ) 

To study the interaction between bounded electrons, the Hamiltonian should be expressed in real space. Eq. ([5|) can 
be transformed into 

Vn=J A e2h2 ^ 2 l%f\^ W (ri) W (r2) ^ (r2) ^ (ri) + ^ (?) 

The k-integral is difficult to calculate for arbitrary ri and r 2 in 3D space. In this work, we study spin-spin interaction 
in a quasi-2D quantum dot sited in x-y plane. Then ri — r 2 sites in x-y plane at z = 0, and k can be integrated out 
from Eq. (JT)) 



0h (k x + ik y ) 2 e^-^-^ _ _27r 2 e 2iarg ( ri - r2 ) 



d V x ^ v > ; = — 3 — , (8) 



where arg(ri — r 2 ) is the angle between ri — r 2 and x axis. So the parallel interaction can be expressed as 

Vb = o ' 2 5 / d 3 rid 3 r 2 - -^4 (ri) ^ (r 2 ) ^ (r 2 ) ^ (n) + h.c. (9) 

32nm z c z eo J |ri-r 2 | 

This equation is valid only when the electrons are confined in a very thin layer. Accordingly, antiparallel interaction 
can be expressed as 



Vex = -7^T^- / d 3 r^(r)^(r)V T W^(r) 



2m 2 c 2 eo 

"T^S— / d'r^r,- - — -^tp\ (ri) ip\ (r 2 ) ^ T (r 2 ) Vj. (ri). (10) 

16m^c^e 7rJ |ri — r 2 | 

The antiparallel interaction conserves the total spin. From this point of view, it is similar to the exchange interaction. 
But there is difference between them. The exchange interaction need the overlap of wave functions, while antiparallel 
interaction need not. 
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III. SPIN-SPIN SCATTERING OF BOUNDED ELECTRONS IN QUANTUM DOTS 

We consider a quasi-2D round quantum dot with radius R. The Hamiltonian is 

H = H + V a - S1 (11) 

where Hq — £ ^ c \, <t c a.<t is the bounded energy and V 8 - a the spin-spin interaction. Because the number of 

electrons in quantum dots is definite, the Coulomb charging energy is a constant and is omitted. We assume the 
confining potential is zero inside the quantum dot and infinite outside it. It is well known that the bounded wave 
function can be expressed as Bessel functions in polar coordinates 

Vw(p,0) = 1 J m e 4 "* (12) 

where m = 0, ±1, • • • is the angular momentum in z direction and Z m ^ n the nth zero of Bessel function J m . The 
bounded energy is 

h 2 

£m ' n = 2^& Z ™' n - (13) 

The spin-spin interaction matrix is calculated 

I r g2iarg(ri-r 2 ) 

(m 3 n 3 | ,m 4 n 4 [\V D \mini ],m 2 n 2 T> = -~r^ \ d 3 rid 3 r 2 — -*-x 

x [ft (r 2 ) (ri) V>i (ri) V 2 (r 2 ) - V> 4 * (ri) V> 3 * fa) V>i (r x ) V> 2 (r 2 )], (14) 

where ^ denotes the wave function tp murH i n Eq. (jT3J) - We work in the atomic units, where h = e = m = 47re = 1. 
By substituting Eq. (fT2"|) into Eq. (Q3]), we obtain 

(m 3 n 3 l,m 4 n 4 l\V D \mini \,m 2 n 2 \) = ^gjfs > ( 15 ) 

where c = 1.370 x 10 2 is the light velocity and a coefficient which depends only upon the quantum numbers m 
and n of initial and final states. 70 can be expressed as 

^ r \ r 2ir 



id = jyyi rFYJ 7V~\ / dpidp2 / d9ld62 

J\m 1 \ + l(^l)-J\m 2 \ + l( Z '2)J\m3\+l{ Z '3)J\m i \+l{^4) Jo Jo 

{pl + pl „ 2piP2CO<ei _ e2))m 1 J - (Z 3Pl )J mi (Z 4P2 )e 

-J m:i (Z3 P 2)J m AZm)e- l{m392+mi9l} l (16) 

where ip = arg(pie 1 — p2 e 2 ) and Zi = Z{m^ ni). Above integral is calculated numerically by Monte Carlo method. 
The 723 is nonzero only when mi,m 2 ,m3 and m 4 satisfy the relation 

7713 + m 4 — mi — m 2 = 2. (17) 

This can be explained by the conservation of total angular momentum. The spin-spin interaction satisfies the com- 
mutation relation [Vd, L z + S z ] = 0. And the bounded states are the eigenstates of L z and S z . So the matrix element 
is not zero only when Eq. (fT7|) is satisfied. We also find that 70 satisfying the selection rule has the magnitude 
of 10 2 ~ 10 3 . For example, j D = 3.823 x 10 3 when the initial and final states are |01 |, -H T> and |01 J., 11 |) 
respectively. 

The second term of antiparallel interaction in Eq. (jTOj) denotes the interaction between two electrons at different 
position, which is proportional to \r\ — r 2 |~ 3 and decreases sharply with the distance of two electrons increasing. The 
size of a quantum dot is usually much larger than the length unit of Bohr radius 0.529A So the second term is usually 
much smaller than the first term and can be neglected. The antiparallel interaction in atomic units is then 

Vex = f d 3 r^ij] (r) Vj (r) ^ (r) ^ (r). (18) 
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The antiparallel interaction matrix is calculated 

(7773773 1,1714114 l\V ex \mim T,m 2 n 2 {) 



d A r- 



-2tt 



^(r)V|(r)Vi (r)^2 (r). 



(19) 



Here we must consider the finite height of the quantum dot. The wave function in circular cylindrical coordinates is 
(choosing the lowest level in z direction [l2j) 



ip(p,9,z) = 



1 



- J 77 



, n P 



a 



) V R 

where a is the height of the quantum dot. Substituting Eq. (I20[) into Eq. ((T9|) . we get 

(7773773 t,m. 4 n 4 UVymmi T,m-2"2 I) = fl.T , 
where 7 ex is a coefficient depending only upon the quantum numbers m and n of initial and final states 

_ Jq dpp,J mi (Z 1 p) J m2 (Z 2 p) J m3 i Z 3P) Jm 4 ( Z iP) 
^|mi| + l (Zl) J\m 2 \ + 1 {^V ^|m 3 | + l C^s) ^|m 4 |+l 

Due to the conservation of angular momentum, the "f ex is not zero only when next selection rule is satisfied 

mi + m 2 — 777,3 — 777,4 = 0. 



(20) 



(21) 



(22) 



(23) 



The 7 ex satisfying the selection rule has the magnitude of 10°. For example, j ex = 0.718 when the initial and final 
states are 1 01 j, —11 |) and | — 11 f, 01 J.) respectively. 



IV. THE EVOLUTION OF TWO PARALLEL SPINS 

We propose a quantum dot with two parallel spins along z direction. In experiments, this can be realized by 
successively injecting two spins into a quantum dot with even number of electrons [l[ @ . Electrons in two highest 
levels are unpaired, while the other levels are occupied by paired spins. So one can assume there are only two 
electrons, occupying two nearby levels. The spins will flip together due to parallel interaction, when the final and 
initial states are degenerate and the transition satisfies the selection rule of Eq. (fTT|) . If the initial and final states 
are non-degenerate, the energy difference is usually much larger than the spin-spin interaction energy. So the non- 
degenerate transition rate is very small and can be ignored. We assume the spins are up at initial time. Arbitrary 
initial state \m\n\ 1,rn 2 n 2 |) has at most three degenerate states | — ttt-itt-i |, 7772772 |), \mini i,—m 2 n 2 J.) and 
I — mini I,— 777 2 772 |). They do not always satisfy the selection rule of Eq. (fT?]) . Only when \rn\\ and \m 2 \ are 
consecutive numbers or one of them is equal to —1, there exists one final state satisfying the selection rule. Below is 
a table showing the quantum numbers m and 77 of several levels sorted by energy. From this table, we can find the 
initial states permitting spin flip. 



level 


1 


2 


3 


4 


5 


6 


7 


8 


9 


m 





1 


-1 


2 


-2 





3 


-3 


1 


n 


1 


1 


1 


1 


1 


2 


1 


1 


2 


ener gy 


5.783 


14.68 


14.68 


26.37 


26.37 


30.47 


40.71 


40.71 


49.22 



For example, the degenerate state of the initial state 1 01 j, —11 j) is 1 01 j, 11 j). In this situation, the initial and final 
states generate the complete bases. Then the Hamiltonian can be written in a matrix 



H = H + V D = 



E V D 
V D E 



where E = {Z^ nl nl + Z^ ll nl )/(2m*R 2 ) is the total bounded energy in atomic units. And 



Vd = <oi i, 11 1 !Vd|01 t, -11 T> 



-Id 



4tt 2 c 2 R 3 



(24) 



(25) 
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is the parallel interaction energy. The evolution matrix evaluates 



U(t) = e- lHt = e~ lEt 
The time-dependent state can be expressed as 



cos Vht —i sin Vot 
-isinVbt cosVoi 



(26) 



|*(t)) = tf(t)|#(0)) = e- lJit cosVbt| TO-ie-^siiiVbtl U>, (27) 

where the quantum number to and n are omitted. The total spin at z direction is oscillating S z (t) = cos2Voi. The 
spins in a quantum dot will precess even without external magnetic field or nuclear spins. The precession frequency 
Vd/tt depends on the size of the quantum dot, being inversely proportional to cube of radius. For example, in a 
quantum dot with radius R — 50nm the precession period is T = tt/Vd — 12.4/is. When it is much shorter than 
the spin relaxation time, the precession can be observed. In GaAs/AlGaAs 2DEG quantum dots, it is impossible to 
observe the precession because the spin relaxation time is only several nanoseconds due to hyperfine interaction [Io| . 
In a quantum dot made of the material with zero nuclear spin (such as Si), the spin relaxation time is much longer. 
Then it maybe possible to observe this precession. 

A weak vertical magnetic field depresses hyperfine interaction and increases the spin relaxation time, but it also 
destroys the precession caused by parallel interaction. Because the degenerate transition becomes the non-degenerate 
transition, whose rate is much smaller, when a magnetic field is added. But if the magnetic field is so strong that the 
Zeeman splitting is comparable with the level spacing, the non-degenerate transitions may become degenerate. And 
the strong field will increase the spin relaxation time. The electron's wave function is more compact in the presence 
of magnetic field, that increases the spin-spin interaction. So the influence of a strong magnetic field to the precession 
is complicated and will be studied in future. 



V. THE EVOLUTION OF TWO ANTIPARALLEL SPINS 



Next we discuss the precession of two antiparallel spins in the antiparallel interaction V ex . We propose a quantum 
dot with two antiparallel spins. At initial time, the two spins occupy two nearby levels |mini "\, mini j). There are at 
most seven degenerate final states, | — mini l\ TO 2"-2 I), \m\n\ tj — m 2"2 I), \ —m\n\ f , — TO2ri2 I), \m\n\ I,m2ri2 f), 
| — mini I, mini T)j |niini !, —min^ X), \ — mini J., —mini f). But the only one permitted by the selection rule of 
Eq. (|23p is |TOini !,TO2n2 |)- This transition is a degenerate transition for arbitrary initial state. The Hamiltonian 
of electrons is 



H = 



E V ex 
V ex E 



(28) 



where V e . 



is the antiparallel energy. The diagonal term of antiparallel interaction is omitted because it is 



c 2 R 2 a 

much smaller than E. The evolution matrix evaluates 



U(t) 



-iEt 



cos V ex t 
-i sin V ex t 



-ismV ex t 
cos V eT t 



The state at time t is expressed as 



e- %Et cos 14**1 Tl> - ie- lEt smV ex t\ |T>- 



(29) 



(30) 



The spins precess synchronously. The precession frequency V ex /2ir depends on the size of the quantum dot, being 
inversely proportional to the volume of the quantum dot. For a quantum dot with the radius of R = 50nm and the 
height of a = 5nm, the precession period is T = 56/is when the initial state is 1 01 f, —11 J.). A magnetic field along z 
direction will not destroy the precession since the initial and final states keep degenerate. 



VI. CONCLUSION 



We have studied spin-spin interaction between two bounded electrons in a quantum dot. The interaction Hamilto- 
nian includes two terms, the parallel interaction and antiparallel interaction, which are shown in Eq. @ and Eq. (|10[) 
respectively. They cause two parallel or antiparallel spins to precess synchronously. If the two spins are parallel, and 
|mi| and |to2| are consecutive numbers or one of them is —1, the total spin oscillates with a period proportional to 
the cube of radius of the dot. A weak magnetic field along z direction will destroy this precession. If the two spins are 
antiparallel and occupy different levels, they always precess synchronously. The period is proportional to the volume 
of the dot. 
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